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<^ Abstract. We study a cosmological model composed of a dark energy fluid interacting with a 

viscous matter fluid in a spatially flat Universe. The matter component represents the baryon and 
dark matter and it is taken into account, through a bulk viscosity, the irreversible process that the 
matter fluid undergoes because of the accelerated expansion of the universe. The bulk viscous 
coefficient is assumed to be proportional to the Hubble parameter. The radiation component is also 
■ taken into account in the model. The model is constrained using the type la supernova observations, 

the shift parameter of the CMB, the acoustic peak of the BAO and the Hubble expansion rate, to 
constrain the values of the barotropic index of dark energy and the bulk viscous coefficient. It is 
found that the bulk viscosity is constrained to be negligible (around zero) from the observations 
Qh, and that the barotropic index for the dark energy to be negative and close to zero too, indicating a 

phantom energy. 
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<N i INTRODUCTION 

> 



In the last years, the type la supernovae (SNe la) observations have given a strong 
evidence of a present accelerated expansion epoch of the Universe (see for instance 
[1, 2, 3] and references therein). 

Several models have been proposed to explain this recent acceleration, one of the most 
successful one is the so-called A Cold Dark Matter (ACDM) that proposes the existence 
of a new kind of component in the Universe called "dark energy" with a behavior of a 
cosmological constant and that constitutes ~ 73% of the total content of matter-energy 
in the Universe today, in addition to a dark matter component filling the Universe in a 
-23% [3]. 

However, this model faces several strong problems, one of them is the huge discrep- 
ancy between its predicted and observed value for the dark energy density (of about 
120 orders of magnitude) [4, 5,6], another one is the so-called the "cosmic coincidence 
problem": the model predicts that we are living in a moment when the matter density in 
the universe is of the same order of magnitude than the dark energy density [7]. 

On the other hand, cosmological models with interacting dark components have been 
studied by several authors, because it is expected that the two dominant components 
(dark energy and matter) interact each other in some way. It has been found that these 
models are promising mechanisms to solve the ACDM problems (see [8, 9] and refer- 
ences therein). 

In addition, it has been known since several years ago before the discovery of the 



present acceleration that a bulk viscous fluid can produce an accelerating cosmology 
(although it was originally proposed in the context of an inflationary period in the early 
universe) [10, 11, 12, 13, 14, 15]. 

So, it is natural to think of the bulk viscous pressure as one of the possible mechanism 
that can accelerate the universe today (see for instance [16, 17, 18, 19, 20, 21]). However, 
this idea faces the problem of that it is necessary to propose a viable mechanism for the 
origin of the bulk viscosity, although in this sense some proposals have been already 
suggested [22, 23]. 

In the present work, following the idea of Kremer et al (201 1) [9] and using the SNe 
la, the shift parameter R of the cosmic microwave background radiation (CMB), the 
baryon acoustic oscillation (BAO) and the Hubble expansion rate H(z) data, we test 
an interacting dark sector model taking into account dissipative process through a bulk 
viscosity in the matter (baryon and dark matter) component, where the interaction term 
is written in terms of the barotropic index of the dark energy fluid. 

In section 1 we present the characteristics of the model and the main equations, in 
section 2 we explain the cosmological probes used to constrain the model and in section 
3 we give our conclusions. 



INTERACTING DARK FLUIDS WITH BULK VISCOSITY 

We study a cosmological model in a spatially flat FRW universe, composed of three 
fluids: radiation, matter and a dark energy fluid components. It is assumed the matter 
component as a pressureless fluid, representing the baryon and dark matter, with a bulk 
viscosity and interacting with the dark energy fluid. 

The Friedmann constraint and the conservation equations can be written as 

? 87TG , 

(pr + Pm + Pde), (1) 

= p r + 4//p r , (2) 

= pm + Pde + 3//(p m + Pde+ J Pm + Aie-3//C), (3) 

where (p r ,p m ,Pd e ) are the densities of the radiation, matter and dark fluid components 
respectively, and (p^Pm.Pde) are their corresponding pressures. The equation (3) arises 
from assuming the interaction between the matter and dark fluid components. The term 
—3H£ corresponds to the bulk viscous pressure of the matter fluid, where £ is the bulk 
viscous coefficient. 

The immediate solution of the conservation equation (2) is 

p I (a)=p l0 /a\ (4) 

where a is the scale factor and the subscript zero labels the present values for the 
densities. 



On the other hand, following the idea of Kremer and Sobreiro [9], the conservation 
equation (3) can be decoupled as 

p m + 3Hy^p m = 0, (5) 

Pde + 3//^ e pde=0, (6) 

where it was defined the effective barotropic indexes and y^ e so that they are related 

as 

r Pm 

with r = p m / Pde corresponds to the ratio between the matter to dark energy densities 
and pi = {ji — \)pi with fi is the usual constant barotropic indexes of the equation of 
state. 

We consider a bulk viscous coefficient £ proportional to the total matter-energy 
density p t = p r + p m + p de , as 

with a dimensionless constant. This parametrization corresponds to a bulk viscosity 
proportional to the expansion rate of the Universe, i.e., to the Hubble parameter [see eq. 

(1)]. 

Following [9] and [24], we assume that the effective barotropic index for the dark 
energy is given as 

^e = 7de + C0. (9) 

So, using (7) and (9), the effective conservation equations (5) and (6) can be rewritten 

as 

Pm + 3//7mPm = 3Hp de Co + 9H 2 ^ (10) 
p de + 3//y de p de = -3//p de Co, (11) 

where it can be identified the interacting term Q = 3//pd e Co- 

Using the expression (9), the solution of the conservation equation (6) becomes 

PdeW=Pde0/a 3(7deH0) . (12) 

On the other hand, with the eqs. (8) and (9) we can express the equation (7) as 

*=*.-£(*.+ig). 

that using the Friedmann constraint (1) we arrive to 

r 

tfn = 7m--^(pr + Pm + 2p de ). (14) 
Pm 



Inserting the eqs. (4) and (12) at (14) we obtain 



T^-r Co (PrO 2 PdeO 



Pm V aA a 3 (rde+^o) 
With this, the eq. (6) for the matter density becomes 



p m + 3H YmPm - 3HH ^ + 2 ^ de0 , - + p m ) = 0. 



(15) 



fl 4 fl 3(y de +Co) 



(16) 



Dividing to (16) by the present critical density p^ rit = 3Hq/(8kG) with //o the Hubble 
constant, and defining the dimensionless parameter densities Cljo = pio/p^, the eq. (16) 
becomes 



d£l m 3 
+ - 



= 0, 



da a 

or in terms of the redshift z with the help of the relation a = 1/(1 +z), 



(17) 



(l+z)^-3 



1 ( 7m -Co)-Co(^(l+z) 4 + 2^ de0 (l+z) 3 ^+^)] =0, (18) 



where it has been defined Cl m = p m /p° rit . The analytical solution of this ordinary differ- 
ential equation (ODE) for Cl m (z) is 



^m(z)=[(l+z)^l2(l+z) 3 ^+2?o) Co( 4_ 3rm + 3 Co)(am0 + ar0 _ 1) _ (19) 

-3(l+z) 4+3 ^Co(7de-rm + 2Co)^ri) + 

(l +z )3r m ((4_3 rm + 3^)( 2 ^ + ( rde _ rm )a m0 ) + 

+ (37de + 37m-8)C(A€)]]/ 

((7de - 7m + 2Co) (4 - 3y m + 3ft)) • 
So, using the solution (19), the Hubble parameter (1) can be written as 

e\z) = a r o(i+z) 4 + a de0 (i+z) 3(rde+?o) +a m (z), (20) 

where E(z) = H{z)/Hq. In the following we will assume y m = 0, i.e., the matter as a 
pressureless fluid. 



COSMOLOGICAL PROBES 



We compare the model with the following cosmological probes that measure the expan- 
sion history of the Universe, to constrain the values of (ft, y de ). 



Type la Supernovae 

We use the type la supernovae (SNe la) of the "Union2" data set (2010) from the 
Supernova Cosmology Project (SCP) composed of 557 SNe la [3]. The luminosity 
distance <1l in a spatially flat Universe is defined as 

a < r L7\ c(l+z) f z dz' 

d L {z, Co, 7de, #o = „ / , r (21) 
H Jo E(z',Co,7de) 

where "c" corresponds to the speed of light in units of km/sec. The theoretical distance 
moduli ji l for the k-th supernova at a distance Zk given by 



Ai f feCo,7de,#o) = 5 log 
So, the x 2 function is defined as 



Mz, Co, 7de, Hp) 
Mpc 



+ 25. (22) 



r,e,«0) - E ( 2 , (23) 

where /i^ is the observed distance moduli of the k-th supernova, with a standard deviation 
of Ok in its measurement, and n = 557. 



Cosmic Microwave Background Radiation 

We use the WMAP 7-years distance priors release shown in table 9 of [25], composed 
of the shift parameter R, the acoustic scale I a and the redshift of decoupling z*. 
The shift parameter R is defined as 



R = H ^^(l+z*)D A (z*), (24) 
c 

where Da is the proper angular diameter distance given by (for a spatially flat Universe) 

c f z dz' 
DA{z)= OTzmJo £(z',Co,7de)- (25) 

With R we can defined a x 2 function as 

v 2 {r v rT,_( R(Co,rde,Ho)-R ohs \ 2 

ZR-CMBlb0,7de,^0) = I ~ I , (26) 

where i? bs = 1-725 is the "observed" value of the shift parameter and Or = 0.018 the 
standard deviation of the measurement (cf. table 9 of [25]). 
The acoustic scale Ia is defined as 

kDa(z*) 

/a = (1+z*J — p—r-, (27) 



where r s (z*) corresponds to the comoving sound horizon at the decoupling epoch of 
photons, z*, given by 



W+z) da 

V3~Jo a 2 H(a)^/\ + (3Q b0 /4Qyo)a' 



(28) 



where we use Q.yo = 2.469 x 10 5 ^ 2 the radiation, and Q b0 = 0.02255/z~ 2 the baryon 
matter component, as reported by Komatsu et al. 2010 [25]. For z* we use the fitting 
formula proposed by Hu and Sugiyama [26] 

z* = 1048 [1 +0.00124(Q b0 /* 2 r a738 ] [l + gi(a m0 h 2 )^} , (29) 

where 

0.0783 (Q h0 h 2 )- - 23 * 0-560 
81 l+39.5(Q b0 /* 2 ) a763 ' 82 l+21.1(ft b0 fc 2 ) 1 - 81 " } 

The X 2 function using the three values (Ia 7 R,z*) is defined as 

3 

XcMB(C0,7de,tf ) = £ (Xi-diXC-^ijiXj-dj), (31) 



where = (Ia,R,z*) are the predicted values by the model and d{ = (Ia = 302.09,7? 



1.725, z* = 1091.3) are the observed ones and C- 1 is the inverse covariance matrix [25] 



2.305 29.698 -1.333 
C" 1 = | 29.698 6825.27 -113.180 ] . (32) 
1.333 113.180 3.414 



Baryon Acoustic Oscillations 



We use the baryon acoustic oscillation (BAO) data from the SDSS 7-years release 
[27]. The distance ratio d z at z = 0.275 is defined as 



«0.275 = 



D v (0.275)' 



(33) 



where zd is the redshift at the baryon drag epoch computed from the following fitting 
formula [28] 



z d = 1291 



(Q. m oh 



2\0.251 



l+0.659(Q m0 /* 2 ) - 828 



l+bi(a m0 h 2 )^ 



b l =03U{Q. m0 h 2 )- QM9 
b 2 = 0.23S(a m0 h 2 f 223 . 



i+o.60i(a m0 h z ) 



2\0.674 



(34) 

(35) 
(36) 



For a flat Universe, Dy(z) is defined as 



D v (z) =c 



dz' 



o H(z') J H(z 



I 1/3 



(37) 



contains the information of the visual distortion of a spherical object due the non 
Euclidianity of the FRW spacetime. 

The <io.275 contains the information of the other two pivots, <io.2 and do.35 usually used 
for other authors, with a precision of 0.04% [27]. 

The X 2 function for BAO is defined as 



*BAo(C0,7de,#0) 



^0.275 — d% b £ 



0.275 



(38) 



where ^275 = 0.139 is the "observed" value and Od — 0.0037 the standard deviation of 
the measurement [27]. 



Hubble expansion rate 

For the Hubble parameter we use 13 available data, 11 comes from the table 2 of 
Stem et al. (2010) [29] and the 2 following data from Gaztanaga et al. 2010 [30]: 
H(z = 0.24) = 79.69 ± 2.32 and H(z = 0.43) = 86.45 ±3.27 km/s/Mpc. For the present 
value of the Hubble parameter we take that reported by Riess et al 201 1 [3 1] H(z = 0) = 
H = 73.8 ± 2.4 km/s/Mpc. The % 2 function is defined as 

xl , M) Ji(^^^, (39 ) 

where H(zi) is the theoretical value predicted by the model and H° bs is the observed 
value. 



Local Second Law of Thermodynamics 

The law of generation of local entropy in a fluid on a FRW space-time can be written 
as [32, 33] 

rv V 5 v = C(v vM v ) 2 = 9// 2 C, (40) 

where T is the temperature and V v s v is the rate of entropy production in a unit volume. 
With this, the second law of the thermodynamics can be written as 

TV V 5 V > 0, (41) 
so, from the expression (40), it simply implies that £ > 0. 
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FIGURE 1. Confidence intervals for (foj/de)- The left panel corresponds to the constraint using the 
SNe la "Union2" data set. The central panel, using the (R,Ia,z*) data from the CMB and the left panel to 
the Hubble parameter H (z) . It was assumed a value of Hq = 73 . 8 km/s/Mpc as suggested by [31]. The best 
estimated values and the % 2 minimum values are shown in table 1. The contours correspond to 63.8%, 
95% and 99% of confidence level. 
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FIGURE 2. Confidence intervals (CI) together for (£oj 7de)- They correspond to the constraints when it 
is used the SNe la dataset, the t/o.275 BAO probe (see eq. [38]), the shift parameter R of the CMB (see eq. 
[26]), the three (R,Ia,z*) data from the CMB (see eq. [31]), the Hubble parameter H(z) (see eq. [39]), and 
the joint "SNe + CMB + BAO + H (z)" (black contours). It was assumed a value of H = 73.8 km/s/Mpc 
as suggested by [31]. The contours correspond to 63.8%, 95% and 99% of confidence level. 



For the present model this inequality becomes (see eq. [8]) 

Co > 0. (42) 



CONCLUSIONS 

It has been studied a cosmological model composed of a bulk viscous matter fluid inter- 
acting with a dark energy fluid. The model is compared with cosmological observations 
to estimate and constrain the values of the bulk viscous coefficient Co proportional to 
the Hubble parameter, and the barotropic index of the dark energy jde- It lS aiS0 use d 



TABLE 1. Best estimated values of the dimensionless coefficients £o and Yde- Figures 1 and 
2 show the confidence intervals. 



Probe 


Co 7de 


■^min ^d.o.f. 


SNe la 

(*,/a,z») CMB 
H(z) 

SNe la + CMB + BAO + H(z) 


-0.102 ±0.022 -0.904 ±0.11 
-0.0025 ±0.001 -0.226 ±0.04 
-0.003 ±0.022 -0.199+j^ 
-0.00035 ±0.0007 -0.363 ±0.028 


590.69 1.06 
1.01 1.01 
8.049 0.731 
635.78 1.06 



the local second law of thermodynamics (LSLT), that states Co > 0, as a criterion for the 
allowed values for Co- 

It is found that using the combined SNe + CMB + BAO + H(z) data sets, the best 
estimated value of Co is negative (implying a violation of the LSLT) and very close 
to zero. The confidence intervals constrain the values of Co to be very around to zero, 
—0.003 < Co < 0.0025, with a 99% of confidence level. We interpret these results as 
an indication that the cosmological data prefer a model with a practically null bulk 
viscosity. Since in the present model the interacting term is proportional to the bulk 
viscosity, this implies also a negligible interaction between the dark components. 

On the other hand, it is found negatives values of y^ e with a 99% of confidence level, 
corresponding to a phantom dark energy. 

It may be an indicative of the phantom energy as a preferred mechanism by the 
cosmological observations (in combination with the LSLT) to explain the accelerated 
expansion of the Universe, instead of the bulk viscous mechanism. 
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